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Abstract 



Boussinesq systems of nonlinear partial differential equations are fundamental equations in 
geophysical fluid dynamics. In this paper, we use asymmetric ideas and moving frames to solve 
the two-dimensional Boussinesq equations with partial viscosity terms studied by Chae {Adv. 
Math. 203 (2006), 497-513) and the three-dimensional stratified rotating Boussinesq equations 
studied by Hsia, Ma and Wang (J. Math. Phys. 48 (2007), no. 6, 06560). We obtain new 
families of explicit exact solutions with multiple parameter functions. Many of them are the 
periodic, quasi-periodic, aperiodic solutions that may have practical significance. By Fourier 
expansion and some of our solutions, one can obtain discontinuous solutions. In addition, Lie 
point symmetries are used to simplify our arguments. 

1 Introduction 

Both the atmospheric and oceanic flows are influenced by the rotation of the earth. In 
fact, the fast rotation and small aspect ratio are two main characteristics of the large 
scale atmospheric and oceanic flows. The small aspect ratio characteristic leads to the 
primitive equations, and the fast rotation leads to the quasi-geostropic equations (cf . [2] , 
[6], [7], [9]). A main objective in climate dynamics and in geophysical fluid dynamics 
is to understand and predict the periodic, quasi-periodic, aperiodic, and fully turbulent 
characteristics of the large scale atmospheric and oceanic flows (e.g., cf. [4], [5]). 
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The Boussinesq system for the incompressible fluid follows in M 2 is 

u t + uu x + vu y — vAu = —p x , v t + uv x + vv y — uAv — 9 = —p y , (1.1) 

9 t + uB x + vO y - kA6 = 0, u x + v y = 0, (1.2) 

where (u, v) is the velocity vector field, p is the scalar pressure, 9 is the scalar temperature, 
v > is the viscosity and k > is the thermal diffusivity. The above system is a simple 
model in atmospheric sciences (e.g., cf. [8]). Chae [1] proved the global regularity, and 
Hou and Li [3] obtained the well-posedness of the above system. 

Aonther slightly simplified version of the system of primitive equations is the three- 
dimensional stratified rotating Boussinesq system (e.g., cf. [7], [9]): 



u t + uu x + vuy + wu z — — v = a(Au — p x ), (1.3) 

Ro 

v t + uv x + Wy + wv z + -^-u = a(Av - p y ), (1.4) 

w t + uw x + vw y + ww z — <jRT = a(Aw — p z ), (1.5) 

T t + uT x + vT y + wT z = AT + w, (1.6) 

u x + v y + w z — 0, (1.7) 



where (u, v, w) is the velocity vector filed, T is the temperature function, p is the pressure 
function, a is the Prandtle number, R is the thermal Rayleigh number and Rq is the Rossby 
number. Moreover, the vector (l/Ro)(— v, u, 0) represents the Coriolis force and the term 
w in (1.6) is derived using stratification. So the above equations are the extensions 
of Navier-Stokes equations by adding the Coriolis force and the stratified temperature 
equation. Due to the Coriolis force, the two-dimensional system (1.1) and (1.2) is not a 
special case of the above three-dimensional system. Hsia, Ma and Wang [4] studied the 
bifurcation and periodic solutions of the above system (1.3)- (1.7). 

In [10], we used the stable range of nonlinear term to solve the equation of nonsta- 
tionary transonic gas flow. Moreover, we [11] solved the three-dimensional Navior-Stokes 
equations by asymmetric techniques and moving frames. Based on the algebraic charac- 
teristics of the above equations, we use in this paper asymmetric ideas and moving frames 
to solve the above two Boussinesq systems of partial differential equations. New families 
of explicit exact solutions with multiple parameter functions are obtained. Many of them 
are the periodic, quasi-periodic, aperiodic solutions that may have practical significance. 
Using Fourier expansion and some of our solutions, one can obtain discontinuous solutions. 
The symmetry transformations for these equations are used to simplify our arguments. 
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For convenience, we always assume that all the involved partial derivatives of related 
functions always exist and we can change orders of taking partial derivatives. The pa- 
rameter functions are so chosen that the involved expressions make sense. We also use 
prime ' to denote the derivative of any one-variable function. 

Observe that the two-dimensional Boussinesq system (1.1) and (1.2) is invariant under 



the action of the following symmetry transformation: 

T(u) = a _1 eiw(a 2 (t + b), ae\{x + a), ae 2 (y + (3)) — a', (1.8) 

T (v) = a- l e 2 v(a 2 {t + b), a(x + a), a(y + (3)) - (3\ (1.9) 

Tip) = a~ 2 p(a 2 (t + b),ae 1 (x + a),ae 2 (y + [3)) + a"x + (3"y + 7, (1.10) 

T{6) = a- 3 e 2 9(a 2 ( y t + b),ae 1 ( y x + a),ae 2 ( y y + P)), (1.11) 



where a, b G K with a ^ 0, €±,e 2 G {1, —1} and a, (3, 7 are arbitrary functions of t. The 
above transformation transforms a solution of the equation (1.1) and (1.2) into another 
solution with additional three parameter functions. 

Denote x = (x, y) . The three-dimensional stratified rotating Boussinesq system is 
invariant under the following transformations: 

%[{u,v,w)\ = ((u(t + b,xA,ez),v(t + b,xA),ez)A,ew), (1.12) 

T 1 (p)=p(t + b,xA,ez), T 1 (T)=T(t + b,xA,ez); (1.13) 

T 2 (u) — u(t, x + a, y + P, ^ + 7) - a', %.{v) = v(t, x + a, y + (3, z + 7) - f3', (1.14) 

T 2 (w) — w(t,x + ot,y + P, z + 7) - 7', %{T) — T(t,x + ot,y + f3, z + 7) - 7, (1.15) 

T 2 (p) = p(t, x + a,y + (3, z + 7) + a~\a" 'x + (3"y + 7" ' z) - R^z + //; (1.16) 

where e = ±1, b G K, A G 0(2, R), and a, /3, 7, // are arbitrary functions of t. The above 
transformations transform a solution of the equation (1.3)-(1.7) into another solution. In 
particular, applying the transformation T 2 to any solution in this paper yields another 
solution with extra four parameter functions. 

To simplify problems, we always solve the Boussinesq systems modulo the above cor- 
responding symmetry transformations, which is an idea that geometers and topologists 
often use. 

The paper is organized as follows. In Section 2, we solve the two-dimensional Boussi- 
nesq equations (1.1)-(1.2) and obtain four families of explicit exact solutions. In Section 
3, we present an approach with u,v,w,T linear in x, y to the equations (1.3)-(1.7), and 
obtain two families of explicit exact solutions. Assuming u z = v z = w zz = T zz = in 
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Section 4, we find another two families of explicit exact solutions of the equations (1.3)- 
(1.7). In Section 5, we obtain a family of explicit exact solutions of (1.3)-(1.7) that are 
independent of x. The status can be changed by applying the transformation in (1.12) 
and (1.13) to them. 

2 Solutions of the 2D Boussinesq Equations 

In this section, we solve the two-dimensional Boussinesq equations (1.1)-(1.2) by an asym- 
metric method and by an moving frame. 

According to the second equation in (1.2), we take the potential form: 

u = ty, v = -£ x (2.1) 
for some functions £ in t, x, y. Then the two-dimensional Boussinesq equations become 

Cyt + CyCxy ~ ^x^yy ~~ ^A^y = ~Vxi £xt + CyCxx ~ CxCxy ~ vA.^ + & = Vyi (2-2) 

Ot + iyOx ~ ixOy - = 0. (2.3) 
By our assumption p xy = p yx , the compatible condition of the equations in (2.2) is 

(AO* + £v(Aflx - U*Z)v ~ ^ + 0x = 0. (2.4) 

Now we first solve the system (2.3) and (2.4). 
Our asymmetric approach is to assume 

9 = e(t,y), Z = <i>{t,y)+xil>{t,y) (2.5) 

for some functions e,<f> and ip in t,y. Then (2.3) becomes 

E t - fay - KEyy = 0. (2.6) 

Moreover, (2.4) becomes 

(Pyy t + Xljjyyt + ((j)y + Xl/jy^yy ~ l/j ((f)yyy + X^yyy) ~ h ' ((f) 'yyyy + X^yyyy) = 0, (2.7) 

equivalently, 

<Pyyt + (fry^Pyy ~ ^<Pyyy ~ V<\>yyyy = 0, (2.8) 
4>yyt + i>yi>yy ~ Wyyy ~ ^yyyy = 0. (2.9) 

The above two equations are equivalent to: 

<t>yt + 0jA - ^<Pyy ~ V<t>yyy = "l, (2-10) 
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i>yt + 4>y - Wyy - v4>yyy = «2 (2.11) 

for some functions a.\ and a 2 of t to be determined. 

Let c be a fixed real constant and let 7 be a fixed function of t. We define 

Ci(s) = 2 ' 771 = 2 ' ( } 

Co(s) = sin 7s, r? (s) = cos 7s. (2-13) 

Then 

+ (-l) r C 2 (^) = c r (2.14) 

and 

9.(Cr(s)) = IVr(s), d s ( Vr (s)) = -(-iy^ r ( S ), (2.15) 

where we treat 0° = 1 when c = r = 0. First we assume 

*P = Piy + foUy) (2.16) 

for some functions (3i and (5 2 of t, where r = 0, 1. Then (2.11) becomes 
(3[ + c r (3 2 2l 2 + (3f + [{(3 2l )' + (-l)V/? 27 3 + 2/? 1 a7 ]i ?r (y) 

+(-i)^ 27 (A 7 - 7')yC,(y) = « 2 , (2.17) 

which is implied by the following equations: 

(3' 1 + c r (3 2 2l 2 + (3l = a 2 , /3 l7 -y = 0, (2.18) 

(/3 27 )' + (-l)V/5 27 3 + 2/3^27 = 0. (2.19) 
For convenience, we assume 7 = \fa' for some function a of t Thus we have 

A = 1 = ^ ? /3 2 = X " , 61 GR. (2.20) 



7 2a" " z ^ 



a 

To solve (2.10), we assume 

= AWrfo) (2.21) 

for some function (3 3 , modulo the transformation in (1.8)-(1.11). Now (2.10) becomes 

[(-l) r ((A»7) ; + A/5s7) + ^3l 3 }Cr(y) = (2.22) 

which is implied by 

(-l) r ((/3 3 7)' + /3i&7) + ^37 3 = 0. (2.23) 

Thus 



where b 2 is a real constant. 

In order to solve (2.6), we assume 

e = be lir,r{y \ 

where b is a real constant and 7 i is a function of t. Then (2.6) is implied by 

iMv) + (-l) r /3 2 77iC%) + K7 2 7i((-l) r ^(?/) - 7iC%)) = 0, 
which is implied by 

i x + (-I)Wti = 0, (-l) r /3 2 - k 77i = 0. 
Then the first equation implies 

7l = b 3 e~ ( - 1)rKa 

for some constant 63. By the second equations in (2.20) and (2.27), we have: 



(-1)' 

For convenience, we take 



VOL 



(-iy uie ^— = b^'e-^ r ^. 
VW 3 



61 = (-l) r ^/c63, &o e R. 

Then (2.29) is implied by 

a / e (-l)>-K)a/2 = ^ 

If u = k, then we have ct = + c . Modulo the transformation in (1.8)-(1.11), 
b = 1 and Co = 0, that is, a — t. When v 7^ k, we similarly take 60 = 1 an d 



a = 



\n[(-iy(is-K)t/2 + c ], c G 



V — K 

Suppose v — k. Then 7 = 1 and 
Moreover, 

= 6ex P (6 3 e-(- 1 ) r ^ r ( J /)), 
£ = ^e-t-^VG/) + (-1)^3^-^^^^) 
by (2.5). According to (2.1), 

u = i y = (-lyi-b.e'^X^+hue-^^riy)], 

v = = -(-lr^e-W^y). 



Note 

ut + uu x + vuy - vAu = blu 2 c r e- ( - 1)r2ut x, (2.38) 
v t + uv x + Wy — vAv — 6 = vv y — b ex^ip^e^^ 1 ^ ut rj r {y)) . (2.39) 
By (1.1), we have 

p = bj eMhe-^'vriy^dy - ^z/V^^^cV + Civ)) (2-40) 

modulo the transformation in (1.8)-(1.11). 
Consider the case v ^ k. Then 

7 = v^7= 1 (2.41) 

y/(-lY(u-K)t/2 + Co 

by (2.32). Moreover, 

= b 2 [(-l) r (is - K)t/2 + c^-^riy) (2.42) 
by (2.21) and (2.24). Furthermore, 

* = 4[(-0>^7/2 i, + c„] + <- 1 >"M(- 1 » " ^ + ^ H '- )+3l2 Uv) (2.43) 
by (2.16), (2.20) and (2.30). According to (2.25), (2.28) and (2.32), 

q = foe M(-i) ? >-KW2+co] 2K/(K - , V(?^ (2.44) 

Similarly, we have 



u t + uu x + vu y - uAu = blc r K 2 [(-l) r (u- K)t/2 + c } 4 ' //{K - u)+2 x 

3(i/ - nfx 
lQ[(-l) r (u-K)t/2 + c } 



N c „ iw.. ..u/o I , 12 » ( 2 - 45 ) 



f t + uv x + Wy — vA — 9 = —tp t + ipipy + wpyy — 6 
= _ te M(-i)->-«)</2+«»r'l»-V(») + l hK{K _ „)[(-!)■•(„ - K )t/2 + c ] 2 "/("-"> +1 / 2 C(!/) 

+ 16[(-l)'(p-*)t/2 4-c ] 2 + ~ ^ + ^"^^to)- P-«> 

According (1.1), we have 

p = b J e b ^ r ^-^ 2+c ^ K,(K ^ 

+5(-l)'J 3 «(« - *)[(-!)> - K )i/2 + c ] 2 "'("-"> +1 , r (! / ) (2.47) 



modulo the transformation in (1.8)-(1. 11). 



Theorem 2.1. Let b, 62, &3, c, cq € K and let r = 0, 1. If v — k, we have the solution 
(2.34), (2.36), (2.37) and (2.40) of the two-dimensional Boussinesq equations (1.1)-(1.2), 
where ( r (y) and i] r (y) are defined in (2.12)-(2.13) with 7 = 1. When u 7^ k, we have the 
following solutions of the two-dimensional Boussinesq equations (1.1)-(1.2): 

u = 4[( J^%^ + Co] + (-irM(-D> - km + Co r^W(y) 

-(-iyb 2 [(-lY(u - K)t/2 + c^l^ +1 /\ r {y), (2.48) 

= ( -iy(y - K) y — _ 3K [(-iy(u - K)t/2 + c ] 2u/(K - u)+3/2 (r(y), (2.49) 

4[(-l) r (z/ - «)t/2 + c ] V ; LV ;V ;/ J V ; V ; 

9 is given in (2.44) and p is given in (2.47), where ( r (y) and r) r (y) are defined in (2.12)- 
(2.13) with 7 = [(-l) r \v - re)f/2 + c ]" 1/2 . 

Observe that 

ip = Guy' 1 (2.50) 
is another solution of (2.11). In order to solve (2.10), we assume 



oo 



= ^7^ (2.51) 

i=i 

modulo the transformation in (1.8)-(1.11), where 7$ are functions of t to be determined. 
Now (2.10) becomes 

00 

-Qu liy - 2 - 18Z/722T 1 + YlH ~ < % + 2)(? + 3)(i + ^H^y'- 1 = «i, (2.52) 

i=i 

equivalently, 

7i = 72 = 0, a- = -60z/73, (2.53) 
1^ - v{% + 2){i + 3)(i + 4) 7 , +2 = 0, i > 1. (2.54) 

72i+2 = 2t7 2i _ g i > x / 2 .55) 

121+2 z/(2i + 2)(2i + 3)(2i + 4) - ' v ; 

_ (2s + l)72i + i = 3607? } 

721+3 v{2% + 3)(2i + 4)(2i + 5) ^(2i + 2)(2i + 5)!' *~ ' 1 ' 

Hence 

= 360 V — y - rr, (2.57) 

* i/»(2i + 3)(2i + 5)!' v ; 

where ct is an arbitrary function of t such that the series converges, say, a polynomial in 
t. 
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Thus 



To solve (2.6), we also assume 

oo 

£ = J>?A (2.58) 

i=0 

where Pi are functions of t. Then (2.6) becomes 

oo 

Gvfcy- 1 + + (i + 2) (61/ - (i + l)^)A +2 ]^ = 0, (2.58) 

i=0 

that is, Pi — and 

P' i -(i + 2)(6v + (i + l)K)p i+2 = 0, i>0- (2.59) 

Hence 

^ 8^y 2i 

e = p + y : — , (2.60) 

where P is an arbitrary function of t such that the series converges, say, a polynomial in 
t. In this case, 

u t + uu x + vUy — uAu = — 60ua, (2-61) 

v t + uv x + vVy — vA — 9 = -36 ^V 3 - - £ 9n n i (R I fo 7TT (2 ' 62) 

~? 2*l!Ilr=l( 6l/ + ( 2r - 1 )«) 

According (1.1), we have 

p = 60uax - 18u 2 y~ 2 + py+) — — (2.63) 

fe2*i!(2< + l)IT r=1 (6i/ + (2r-l)«) 

modulo the transformation in (1.8)- (1 . 1 1) . 

Theorem 2.2. We have the following solutions of the two-dimensional Boussinesq 
equations (1.1)-(1.2): 



™ a {i) 2i+2 



i=0 



# zs given in (2.60) and p is given in (2.63), where a and P are arbitrary functions of t 
such that the related series converge, say, polynomials in t. 

Let 7 be a function of t. Denote the moving frame 

w = x cos 7 + y sin 7, w = y cos 7 — xsin7. (2.65) 

Then 

d t {vj) = 7'ro, d t (w) = —^'w. (2.66) 



Moreover, 

= cos 7 d x + sin 7 d y , d& = — sin 7 d x + cos 7 (2.67) 

In particular, 

A = dl + d 2 = d% + d%. (2.68) 

We assume 

i = 0(t, w) - ^-{x 2 + y 2 ), 6 = ^(t, w), (2.69) 
where and ^ are functions in t, w. Then (2.3) becomes 

A - «fe = (2.70) 

and (2.4) becomes 

-27" + <f> t && - v<f>tt,&Tv + i>™ cos 7 = 0. (2.71) 
Modulo the transformation in (1.8)-(1.11), the above equation is equivalent to 

-li'w + <p t{b - v<f>tt& + if? cos 7 = 0. (2.72) 

Assume v — k. We take the following solution of (2.70): 

m 

i> = Y^ a ^e a ? Ktcos2bl+a ^ cosbl sm(a 2 Kt sin 2b { + caw sin k + k + a) (2.73) 
i=i 

where a^, &j, q, di are real numbers. Moreover, (2.72) is equivalent to solving the following 
equation: 

m 

Ivi - i'w 2 + <f> t - + [J2 d l e a > tcos2h+a ^ cosh * 

i=i 

x sin(a 2 nt sin 26j + aiw sin hi + q)] cos 7 = (2.74) 
by (2.1). Thus we have the following solution of (2.74): 

m „ 
= _[^ d . e a?*ctco S 26 i +a i *cos6 i ^(^^2^ + SU1 6; + C,)] / COS 7 dt 

1=1 

+ 7 'w 2 + ^ J se a^tcos2S s +a s ^cosS s sin (a2^ sin 26 s + d s wsmb s + c s ), (2.75) 



s=l 



where d s ,b s ,c s , d s are real numbers. 

Suppose v 7^ k. To make (2.72) solvable, we choose the following solution of (2.70): 

m 

= s ^a l d i e a > t+a ^ . (2.76) 
i=i 



10 



Now (2.72) is equivalent to solving the following equation: 

m 

vi - i'w 2 + <f) t - v^tt + d ie a * Kt+ai ™ cos 7 = (2.77) 

i=i 

by (2.1). We obtain the following solution of (2.77): 

n 

= 7 '^ 2 + ^ d s e^ Ktcos2k+ ^ cosk sm(a 2 s Kt sin 2b s + a s w sin b s + c s ) 

s=l 

m „ 

- die a > t+a ^ / e a ?( K -^* cos 7 dt. (2.78) 



i=l 

Note 

/. .. . 1 



u = ro sin7 - j'y, i; = 7'x - ro cos 7. (2.79) 
By (2.72), 

■u t + uu x + vu y — uAu 
= {4>™t - Ww ro ) sin 7 + 27'0 ro cos 7 - ^' 2 x - 7"?/ 
= {2^''w — ip cos 7) sin 7 + 27'0 ro cos 7 — 7 /2 :r — y'y, 

= 7''(a;sin27 — y cos 27) + (27'0 ro — ■0sin7) cos 7 — 7' 2 x, (2.80) 

i>i + -Ufa; + vv y — uAv — 9 
= ( Ww ro - <t>wt) cos 7 + 27> ro sin 7 - 7 /2 y + 7"x - ip 
= (ip cos 7 — 2y'w) cos 7 + 27' sin 7 — 7' 2 y + 7 / 'x — ip 

= — 7''(a;cos27 + ysin27) + (27'0 ro — ^sin7) sin7 — 7' 2 y. (2-81) 

According to (1.1), 

l' 2 ~ l'' s i n 27 o 7 /2 + 7'' sin 27 2 ./ f . 

p— x H y +7 xy cos 27+ / -^cfaj sin 7 — 27 <p (2.82) 

modulo the transformation in (1.8)-(1. 11). 

Theorem 2.3. Let 7 be any function of t and denote w = a; cos 7 + y sin 7. Take 

{a,, 6j, Q, dj, d s , 6 S , c s , d s | i = 1, m; s = 1, n} C R. (2.83) 

If v = k, we have the following solutions of the two-dimensional Boussinesq equations 
(1.1)-(1.2): 



= -iy + sin7{27'ro + ^ a s d s e% Ktms2k+iis ™ cosk sin(a 2 /tt sin 2b s + a s wsmb s + b s + c s ) 

s=l 

m „ 

-[^a l d i e a ? Ktcos2bl+a ^ cosbl sin(a 2 /ttsin26 J + 6 l + a l rosin6 J + c l )] / cos 7 dt}, (2.84) 
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v = jx- cos7{2 7 'ro + ^a s 4e a ' Kicos2£s+a ^ cosSs sm(a 2 s Kt sin 2b s + a s iusmb s + b s + c s ) 

s=l 

m „ 

J2a i d i e a l Ktcos2bi+ai ™ cosbi sm(a 2 Ktsm2bi + turn nub, + h + a)} / cos 7 dt}, (2.85) 



i=i 

6 = ip m (2. 73), and 



cos 7 )[^ d t e a > tcos2b > +a ^ cosb > sm(a 2 Ktsm2bi + a t zusmh + a)] 
i=i 

7 /2 -7 // sin27 2 7 /2 + 7"sin27 2 , , 2 _ 2 

H x H y +7 it/ cos 27 — 27 w 

2 2 



n 

A2„ 



-27'^4e a ' Ktcos26s+a ^ cos6s sin(a 2 /ttsin26 s + a s rosin6 s + c s ). (2.86) 

s=l 

When u ^ k, we have the following solutions of the two-dimensional Boussinesq equa- 
tions (1.1)-(1.2): 



n 

u = {^a s d s e^ Ktcos2bs+ ~ a ^ cosb ^ 
s=i 

m „ 

+2iw - J2 aAe a > t+a ^ I e a ^ K -^ 1 cos 7 dt} sin 7 - 7 'y, 
i=i ^ 

n 

v = -{^MV' K * cos2Ss+as " cosis sin(a>tsin26 s + a s rosin6 s + 6 s + c s ) 

s=l 

+2iw - J2 a i d i e a * vt+ai * I e a ? (K - iy)t cos7rft}cos7 + 7'x, 
i=i ^ 

= ip in (2. 76), and 

7 /2 - 7" sin 27 2 7' 2 + 7"sin27 2 , , 2 _ 2 

p = x H y + 7 rcy cos 27 — 27 

n 

-2i^d s e^ Ktcos2bs+!l ^ cosk sin(a 2 fi:tsin26 s + a s wsin b s + c s ) 

s=l 

m » 

+ ^rf ie a ^ t+a ^(27' + sin7) e a ^ K - u)t cos -fdt). (2.89) 



(2.87) 



(2.88) 



Remark 2.4. By Fourier expansion, we can use the above solution to obtain the one 
depending on two piecewise continuous functions of w. 

3 Asymmetric Approach I to the 3D Equations 

Starting from this section, we use asymmetric approaches developed in [11] to solve the 
stratified rotating Boussinesq equations (1.3)-(1.7). 

12 



For convenience of computation, we denote 

$ 1 = u t + uu x + vu y + wu z - -^-v - a{u xx + u yy + u zz ), (3.1) 

<& 2 = vt + uv x + vv y + wv z + -^-u - a(v xx + Vyy + v z z ) , (3.2) 

$3 — w t + uw x + vWy + ww z — aRT — a(w xx + w yy + w zz ). (3.3) 
Then the equations (1.3)-(1.5) become 

$i + crp x = 0, $ 2 + ap y = 0, $ 3 + ap z = 0. (3.4) 

Our strategy is to solve the following compatibility conditions: 

= d x (Q 2 ), = d x (Q 3 ), <9 2 ($ 2 ) = ^($3). (3.5) 

First we assume 

u = <j) z (t, z)x + q(t, z)y + fi(t, z), v = r(t, z)x + ip z (t, z)y + e(t, z), (3.6) 

w = -<p{t,z) -i/>(t,z), T = &(t, z) + z, (3.7) 
where 0, 1?, q, fi, r, and £ are functions of t, z to be determined. Then 

$1 = <ptzx + q t y + l^t + (t)z((pzx + <,y + ii) + (q-l/RQ)(Tx + i) z y + e) 
-(0 + i)){<\) zz x + fey + /i 2 ) - a((p zzz x + q zz y + fx zz ) 

= [<t>tz + (fi + r (<? - V^o) - 4>zz{4> + - o-0^J^ 

+Ht + n<p z + (<r - 1/ R )s - Hz{<t> + il>) ~ p^zz, (3.8) 

$2 = T t x + ^ tz y + e t + ^ z (rx + tlj z y + e) + (T + l/R )((f) z x + qy + /j,) 
-(0 + i/))(t z x + ij zz y + e z ) - ct(t zz x + ^ zzz y + e zz ) 

= [iptz + 4> 2 z + ?(r + 1/ i? ) - (0 + i>)^zz - V^zzz]y 
+[r t + rip z + (r + 1/R )(f) z - ((f) + i))r z - ot zz \x 

+s t + eip z + (r + l/Rojii - (0 + ip)e z - ae zz , (3.9) 

$3 = ~4>t ~ A + (0 + ^)(0 2 + ^) - <7i2(tf + 2) + <7(0 Z2 + 4> zz ). (3.10) 
Thus (3.5) is equivalent to the following system of partial differential equations: 

<ptz + <Pl + r(q - 1/Ro) - 4> zz {4> + V) - o<i>zzz = <*i, (3-11) 
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<h + S(f>z + ^z^ - l/Ro) - <^(0 + VO ~ VSzz = «, (3.12) 

Ik + v4z + (<?-!/ ^o)^ - Vz{<t> + ~ °^zz = a 2 , (3.13) 

i> tz + ^ + ?(r + 1/Rq) - (0 + - = (3.14) 

T% + "H^ + (t + l/Ro)4> z - (0 + VO 7 "* - ffr 22 = a, (3.15) 

e t + + (t + 1/-Ro)a* - (0 + - oe st = fa (3.16) 

for some a, a±, a 2 , Pi, 02 are functions of t. 

Let 0^6 and c be fixed real constants. Recall the notions in (2.12) and (2.13) with 
7 = 6. We assume 

= rSCr </> = & -1 (72C^) + l3Vr(z)), (3.17) 

? = 74(72^(2) - (-l) r 73Cr(^)), r = 7571^(2), 7475 = 1, (3.18) 
where 7$ are functions of t to be determined. Moreover, (3.11) becomes 

(Vi + (-l) r &V 7 i - 7175/^0)^(2) + (71 + 72)7ic r = «i, (3.19) 
which is implied by 

<*i = (71 + 7 2 )7ic r , (3.20) 

7 ; + (-l) r 6V 7 i - 7175/^0 = 0. (3.21) 
On the other hand, (3.15) becomes 

[(7i7 5 )' + 7i/^o + (-l) r b 2 a^ 5 }r] r + 7l75 ( 7l + 72 )c r = a, (3.22) 

which gives 

a = 7i75(7i + 72)c r , (3.23) 

(7i7 5 )' + (-l) r 6V 7l75 + 7l /i2o = 0. (3.24) 
Solving (3.21) and (3.24) for 71 and 7175, we get 

7i = & 1 e-<- 1 > r " rt sin A, 7175 = M-^^ cos A (3.25) 

where b\ is a real constant. In particular, we take 

75 = cot (3.26) 

Observe that (3.12) becomes 

[(7274)' + (-l) r 6V 727 4 - 72/^^(2) + 74(7i72 + 72 + (-l) r 7 3 2 )c r 

-(-l) r [(7374) / + (-l) r 6V 72 74 - 73/i2o]Cr(^) = « (3.27) 
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and (3.14) becomes 

[y 2 + (-l)'&Vy 2 + Wa/RMz) + (7i72 + ll + (-l) r 7l)c r 

-(-l) r [73 + (-1)^V 7 3 + 7374/i? ]Cr(^) = Pi, 

equivalently, 

a = 74(7172 + 72 + {- l YllVi 
A = (7i72 + 72 + (-l) r 7 3 V, 
(7274)' + (-l) r fr 2 C7274 - 72/-R0 = 0, 

72 + (-l) r &V 72 + 7274/^0 = 0, 
(7374)' + (-l) r 6 2 o"7274 - 73/-R0 = 0, 

7 3 + (-l) r 6V 73 + 7374/^0 = 0. 
Solving (3.31)-(3.34) under the assumption 7475 = 1, we obtain 

7274 = & 2 e-<- 1 > r * rt sin ' 72 = b 2 e'^ h ^ cos ± 

Kq Kq 

7374 = b 3 e-(-^ S i n * 73 = b^' 1 ^ cos ±-. 

K Ko 

In particular, we have: 

* 

74 = tan — . 

According to (3.23) and (3.29), 

7i75(7i + 72)c r = 74(7i72 + 7 2 + (-l) r 7l)c r , 



equivalently 



Thus 



So 



It 9+ 
-26462 cos — + (b 2 2 - b\ + {-!)%) sin — = 0. 



6462 = 0, b 2 2 -bi + (-l) r bl = 0. 
r = 0, 6 2 = 0, 64 = 63 



or 



r = 1, 61 = 0, 6 2 = 63. 
Assume r = and 61 7^ 0. Then 

= b~ 1 bie~ b2(Tt sin 62 sin ^ = 6~ 1 6 1 e~ b2,T * cos 62 cos 

ito 
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q = — b\e b2(7t sin bz sin r = b x e b2(Tt cosbz cos 

Rq Rq 



Moreover, we take /i — e — $ — 0. Then 

$ 1 = 1 2 ( x + 7 5 y) = 6^e~ 262,Ti sin 



. t t 
. x sin — + y cos — 



by (3.8), (3.11)-(3.12), (3.20) and (3.23). Similarly 

t ( 



t 



$ 2 = b\e 2b2(Tt cos — I x sin — + y cos — ) . 

Kq \ Rq Rq 



(3.44) 



(3.45) 



(3.46) 



According to (3.10) 



b-'R^he'^ - b- x b\e- 2h2 ° l cos I fez - — 

/to 



sin i 62 — — ) — Tin : . 



(3.47) 



By (3.4), we have 

Rz 2 b ie ~ b ' 2at 



P 



2 ' b 2 aR 

b 2 e -2**t 



t \ b(e 

Rq 



2 a -2b 2 crt 



+ -tt; — ^— cos [bz — — ) — " x „ — — cos 2 [bz — — 



la 



t 



2ab 2 
I . It 



t 

Rq 



y 2 cos 2 — + x 2 sin 2 — + xy sin — 

Rq Rq Rq 



(3.48) 



modulo the transformation in (1.14)-(1.16). 
Suppose r = 1 and b 2 7^ 0. Then 

= r = /i = e = tf = O, i) = b'%e bz+b2at cos q = b 2 e bz+b2<jt sin 4 

Rq Rq 

Moreover, 



"ET- (3-49) 



= $ 2 = 0, $ 3 = fe-^ai^ 1 e 62+b2(T * sin + & _1 ^e 2 ( te+62<rt ) cos 2 - Raz. 

Rq Rq 

According to (3.4), 



V 



Rz 2 b 2 e bz+b2at . t % e 2(bz+#<Tt) 2 t 



b 2 aRn Sin Rq 



2 U'OJXq JXq 

modulo the transformation (1.14)- (1.1 6). 



Wa— COS Rq 



(3.50) 



(3.51) 



Theorem 3.1. Let b,b 1 ,b 2 G K twt/j 7^ 0. VFe nave £/ie following solutions of the 
three-dimensional stratified rotating Boussinesq equations (1.3)-(1.7): (1) 

u = bie~ b2at (x cosbz — y sinte) sin t> = bie~ b2 at (x cos bz — y sin bz) cos (3.52) 

/to Rq 



w = -b^he-*" cos ( bz - ^- J , T = z 



and p is given m (3.48); (2) 



u = b 2 e bz+b ' 2at y sin u = b 2 e bz+b2(7t y cos - 

/l ito 



(3.53) 



(3.54) 
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w = -b~ 1 b 2 e bz+b2at cos T = z (3.55) 

and p is given in (3.51). 

Next we assume = <^ = -?/> = t = O. Then 

Ht - -^-e - (T// Z2 = a 2 , e t + -^-i/ - oe zz = fo, &t ~ ^zz = 0. (3.56) 

Solving them, we get: 

Theorem 3.2. Let bi, Ci, di, a r ,b r , c r , d r , a s ,b s , c s , d s be real numbers. We have the 
following solutions of the three-dimensional stratified rotating Boussinesq equations (1.3)- 
(1.7): 

m 

u = cos— Y^d z e a ^ tcos2bMcosb ' sm(a 2 atsm2b i + a i zsmb i + c i ) 
t n 

+ sin — V d r e^ at cos 2br+arZ cos br sm(a 2 at sin 2b r + a r z sin b r + c r ) , (3.57) 



v = -sin— J" d t e a > tcos2b ^ zcosh sm(a 2 atsm2bi + a lZ smb t + c t ) 



i=l 



+ cos^- V J r e a " CTicos2Sr+a ^ cos ^sin(a2atsin26 r + a r zsin6 r + c r ), (3.58) 



k 

IV = 0, T = 2; + ^a s d s e a » icos4+5s2COS ^ sin(a^sin26 s + fi s ,2sin6 s + 6 S + c s ), (3.59) 

s=l 



Hz 2 m ' A 

p = — - + d s e & ' tcos2bs+&sZCOsbs sm(a 2 s tsm2b s + S s 2sin6 s + c s ). (3.60) 



s=l 



Remark 3.3. By Fourier expansion, we can use the above solution to obtain the one 
depending on three arbitrary piecewise continuous functions of z. 

4 Asymmetric Approach II to the 3D Equations 

In this section, we solve the stratified rotating Boussinesq equations (1.4)- (1.7) under the 
assumption 

u z = v z = w zz = T zz = 0. (4.1) 
Let 7 be a function of t and we use the moving frame w in (2.65). Assume 

u — f(t, w) sin 7 — j'y, v — — f(t, w) COS7 + j'x, (4.2) 
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According to (4.3), we assume 

w — (f)(t,zu), T = ip{t,w) + z, (4.3) 

for some functions /, and ip in t and w. Using (2.66)-(2.68), we get 

$i = -{i 2 + l'/Ro)x - i'y + ft sin 7 + (2Y + l/R )f cos 7 - af** sin 7 , (4.4) 

$2 = -(l' 2 + l'/Ro)y + i'x - ft cos 7 + (2V + l/i? )/ sin 7 + cos 7, (4.5) 

$ 3 = &-(7&& 1& -<7.R(^ + z). (4.6) 

By (3.5), we have 

-2 7 " + Ut - af*** = 0, (4.7) 

(fit - ofyixvj - crR?p = 0. (4.8) 

Moreover, (1.6) becomes 

A - = 0. (4.9) 

Solving (4.7), we have: 



rn 



f = 2-f'w + a ^ je a ? Ktcos26l+a ^ cosbl sin(a?Kf sin 26* + a { w sin 6, + h + q), (4.10) 
i=i 

where cti,bi,Ci,di are arbitrary real numbers. Moreover, (4.8) and (4.9) yield 

n 

= ^ d r e^ tcos2hr+ ^ cosir sm{a 2 r tsm2bi + a r wsinb r + c r ) + aitt^, (4.11) 



r=l 

fc 

^ = J se a2tcos26.+a.*ooeS. sin (a2 tsin2 5 s + a s wsin6 s + c s ) (4.12) 

s=l 

if cr = 1, and 

n 

= ^J re a^tcos2S r +a^cosS r gin^a^^gg. + — gin ^ + a r ) 

r=l 

+ 7^-^ ^ J se ^tcos26 3 +a^cosS s ^ 2 6 s + fi s T& Sin 6 a + C s ), (4.13) 

°" s=l 

k 

^ = J2ald s e^ tcos2k+ ' a ^ cosk sm(a 2 s tsin2~b s + d s msmb s + 26 s + c s ) (4.14) 

s=l 

when (j^l, where a r ,b r ,c r ,d r ,a s ,b s ,c s , d s are arbitrary real numbers. 
Now 

$1 = (7 /, sin27-7' 2 -77i?o)a;-7"l/cos27 + (27 / + l/ J R )/cos7, (4.15) 
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$2 = -(7 // sm27 + 7 ,2 + 77^o)y-7 / ^cos2 7 + (2 7 ' + l/i? )/sin 7 (4.16) 
and $ 3 = —oRz. According (3.4), we have 

p = - 27 ' + 1 / R ° \iw 2 + d ^ iKt C ° S 2k+ai ™ C ° S k sin(a?«t sin 2b t + a { w sin ^ + q)] 

i=i 

2 (7 /2 + 7V^o)(^ 2 + y 2 )+7 // (l/ 2 -a: 2 )sin27 7 " 
H — 2; H 1 x|/cos27 (4.17) 

2 2(7 O" 

modulo the transformation in (1.14)- (1.1 6). 

Theorem 4.1. Lei ai,bi,Ci,di,a r ,b r ,c r ,d r ,a s ,b s ,c s ,d s be real numbers and let 7 be 
any function of t. Denote w = a; cos 7 + ysin 7 . We have the following solutions of the 
three-dimensional stratified rotating Boussinesq equations (1.3)-(1.7): 



m 

U 



a ^e a ? Ktcos2fel+a ^ cos6 ' sin(a 2 /tt sin 2b, + a { w sin k + k + a) 



i=i 



+2^'w] sin 7 - 7'y, (4.18) 



u = [-^a^e a ? Kicos2bl+a ^ c ^ 
i=i 

+27'^] C0S7 + 7'x, (4.19) 

n 

w = d r e^ tcos2ir+ ^ cosir sin(a 2 t sin 26< + a r w sin 6 r + c r ) 

r=l 

+(jfflj](f s e s ' (cos2l+s ^ cos ^ sin(5 2 tsin26 s + 5 s rosin6 s + c s ), (4.20) 



8=1 



T = 2 + J se a?tcos2b s +a s &cos6 s ^ 3^ + ^ sin ^ + Q (4.21) 

s=l 

if a = \, and 

n 

d r e^ atcos2hr+ ^ cos ~ br sm(a 2 r atsm2bi + a r wsmb r + c r ) 



w = 

r=l 

k 



+ ^^^^ ea '* C ° S26s+as " COS ^ sin (^ tsin2 ^ + a s ^smb s + c s ), (4.22) 

s=l 

k 

T = z + J2a 2 s d s e^ tcos2is+ ^ cosis sin(a 2 tsin26 s + d s wsinb s + 2b s + c 8 ) (4.23) 



s=l 



when d^l. 
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Remark 4.2. By Fourier expansion, we can use the above solution to obtain the one 
depending on three arbitrary piecewise continuous functions of w. 



Next we let a be any fixed function of t and set 

zu = a(x 2 + y 2 ). (4.24) 

We assume 



a' a' 



u = y(f>(t, w) - —x, v = -x(f)(t, w) - —y, (4.25) 
a' 

w = ip(t,zv) + —z, T = &(t,vj) + z (4.26) 
a 

where 0, ip and d are functions in t, zu. Note 

a' 2 + 2aa'' a' ...fx a' 



*' = Ia^ X + 2R^ + v4 " + U " « v ) * " ** " 4<ra ^>— < 427 » 

$ 2 = ^^y - _* - sfc + + -* ) * - yf + W(^W (4.28) 

According to the first equation in (3.5), 



a' 



a 



a' 



+ 2^ = °' ^ 



equivalently, 
for some function (3 of t. Write 



ot' , a'zu 



™4>)t w 4> - 4craG7(G70) roro + — — = a(3' (4.30) 



^ a 2i? « ^ ^ 



Then (4.30) becomes 

t - 4<jcra0 roro = 0. (4.32) 

Suppose 

oo 

4> = J2^\ (4-33) 
i=i 

where 7$ are functions of £ to be determined. Equation (4.32) yields 

( 7i ) t = 4i(i + l)<T07 i+ i. (4.34) 

Hence 



for some function 7 of t. Thus 



7j+l = wttw^y ( } 
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By (4.31), we get 

= ^-2^ + a ^,!(, + l)!(4^- (437) 

Note 



a' 



$3 = V>t + —i> - Ao{w^ w ) w - aR($ + z). (4.38) 
a 

By the last two equations in (3.5), 

a' 

i/>t + —i/>- 4a(z^ ro ) ro -aR& = (4.39) 
a 

modulo the transformation in (1.14)- (1.1 6). On the other hand, (1.6) becomes 

t -4Mg w = O. (4.40) 



Hence 

<? = v gi ro (4.4i) 

modulo the transformation in (1.14)-(1.16), where 0\ is an arbitrary function of t. Sub- 
stituting (4.41) into (4.39), we obtain 



^ (4cr) J ((i + l)!) 2 



i=i 

where # 2 is another arbitrary function of t. 
Now 



(4.42) 



$1 = — x H h —0 - £0 , (4.43) 



4a 2 zu R 

+ 2cm'' a(3'x y 
~4a* y ~^^ + R~ 



$ 2 = — y -J^ + JL^-ytf (4.44) 



by (4.27) and (4.28), and 

$ 3 = (qtV - <r.R)z (4.45) 
by (4.38). According to (3.4), we obtain 

a' 2 + 2aa'' 3 \ . 2 2 . (3' y (Ro&l — 



P = I 4aa 2 + 8^|J (3: + y ) + ~ arctan - + ^ lna(s+y) 
a- 1 /? 2 oR-a^a'R 2 1 A (cT^'fr^+V + S/ 2 ) <+1 

+ ^ 2 -^rl. 



2(x 2 + y 2 ) 2a ((i + l)!) 2 (4 



(J 



00 



(a- 1 ft) i ( 7 )(a-W(7)(«(^ + y 2 )) 



2^ + l)!(j + 1)K« + J + 1)(4<t)^ 



2a f-f i\{i + l)\i{AaY 1 ' ' 
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modulo the transformation in (1.14)-(1.16). By (4.25), (4.26), (4.37), (4.41) and (4.42), 
we have: 



Theorem 4.3 Let a, /3, 7, #i, 62 be any function oft such that the following involved 
series converge. We have the following solutions of the three-dimensional stratified rotating 
Boussinesq equations (1.3)-(1.7): 

X i + y* 2R 2a + V ^ o + 1)!(4<t)< ' 1 ] 

2R 2a V x i + y* +aX ^ + 1)!(4<t)< ' 1 8J 



a «tT (4^) l ((^ + l)!) 

' 0} l V + V + 2/ 2 ) 



4^(( ? + l)!)2 

an<ip is given m (4-46). 

5 Asymmetric Approach III to the 3D Equations 

In this section, we solve (1.3)-(1.7) with v x = w x = T x = 0. 
Let c be a real constant. Set 

w = y cose + z sine. (5.1) 

Suppose 

u = f(t,w), v — <f>(t, va) sine, (5.2) 
w = — (j)(t, w) cose, T = ip(t,zu) + z, (5.3) 
where /, and ip are functions in t and gj. Then 

$1 = ft ~ Vfvaw ~ (5.4) 

$2 = (<pt - o~(j) mzu ) sinc+ -^-/, (5.5) 
$3 = (o-0 roro - 4>t) cose - <ri?(^ + z). (5.6) 



By (3.5), 



Sill c 

/wt - o"/ rororo 5—0*7 = °» ( 5 - 7 ) 

Sill c 

&t - o-0 roro ) ro + -5— /w + oRtyw cos c = 0. (5.8) 
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Modulo (1.14)- (1.1 6), we have 



sin c 

f t - - —<f> = 0, (5.9) 

Ko 

Sill c 

(fit - o<\>vjvj + -5— / + crRtpcosc = 0. (5.10) 



Denote 



1 I \ oin tsinc t sin c I 1 j 



(5.11) 



Then (5.9) and (5.10) become 



^ £ sin c 

ft- cff^ - aRificosc sm— — = 0, (5.12) 

a £ Sill c 

t — <j(fi muj + cri?^ cos c cos — - — = 0. (5.13) 

Ro 

On the other hand, (1.6) becomes 

V> t -Vw = 0. (5.14) 
Assume a = 1. We have the following solution: 

m 

if, = ^a^ je a ? tcos2fe ' +airocosfel sin(a^sin2fe i + a^sin^ + b> + q), (5.15) 
i=i 



f = -RR COt C COS V a i d i e o ' t00s26i+OiW00B6i sin ( a 2 t gin 2 ft. + sin ft . + ft. + Q ) 

n 

+ ^a r d r e^ tcos2K+ ~ a ^ cosK sm(a 2 r tsm2bi + a r wsmb r + b r + c r ), (5.16) 

r=l 



m 

<f> = -RR COt C Sin V a ^. e «?^os26 l +a iro cosb l ^ 2ft . + ^ ^ ft . + ft . + Q ) 

+ J2^d s e h2jcos2k+ ' aamcosk sm(a 2 s tsm2~b s + a s wsin~b s + ~b s + c s ), (5.17) 

s=l 

where aj, q, q, a r , 6 r , c r , d r ,a s ,b s , c s , d s are arbitrary real numbers. According to (5.11), 



m 



f = -RR cotc cos^^ *y a i d i e e t ta » 2bi+aiW ' xebi sin(a 2 t sin 2b, t + a, t w sin 6, + h + q) 



r sin c 



+ cos — — V a r rf r e a ' tcos2S '' +a '- rocos£ '- sin(a 2 t sin 26< + a r w sin 6 r + 6 r + c r ) 
/t 



r=l 



. fc 

+ S in-^ Y^ Ss J se a^cos26 s +a sro cos6 s ^(^^3^ + ^ gin ft s + ft s + (5.18) 
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. t sin C ^ A 3 . : 
<p — — sm — — > a r d r e 7 



t cos 26 r +a r ra cos fe r 



o 



sin(a 2 t sin 2bi + a r w sin b r + b r + c r ) 



r=l 



t sin c v-^ _ 7 
+ cos — - — y, a sd s 



J a,tcos2fe s +a s rocosb s ■ ( ~2 



sin(a s i sin 2b s + a s zusinb s + b s + c s ). (5.19) 



o 



s=l 



Suppose (7^1. We take the following solution of (5.11)-(5.14): 

m 



(5.20) 



i=i 



./ 



m cose 



af (1 - cr) sin ^ - i? 1 sine cos ^ 



aKl-<r) 2 + i^ 2 sin 2 c 
^a r rf r e a ' atcos2S '' +ai - rocos6 '- sin(a 2 crtsin26 i + a r wsmb r + b r + c r ), 



(5.21) 



r=l 



= fT J R^a i ^e a ?* +ai 
i=i 
k 



cose 



af (cr - 1) cos - i? sin c sin 



t sin c 



at cos 26 s +a s ro cos 6 S 



af(l-a) 2 + R 2 sin 2 c 
sin(a 2 crt sin 26 s + a s cu sin b s + b s + c SJ . 



(5.22) 



s=l 



where ai,bi,Ci,a r ,b r ,c r ,d r ,a s ,b s ,c s , d s are arbitrary real numbers. According to (5.11), 
t sin c 



t bill C ^— \ , ' fi 2 

/ = cos — — > a r d r e 1 
Rr — ' 



crt cos 2fc r +a r ro cos ? 



sin(a r crt sin 26j + a r trj sin b r + b r + c r ) 



r=l 



+ sin a^e^ tcos2 ^ +(isrocos ^ sin^crt sin26 s + 5 s wsin6 s + b s + c s 



Ro 



s=l 



aidie°Z t+a ^ sin 2c 



^ 2^(^(1 - cr) 2 + # " 2 sin 2 c) ' 



(5.23) 







t sine 



— sm • 



Ro 



^^^^ Qj^d^G 



r<rt cos 2b r +a r ro cos b r „: T1 /'ft2 



sin(a r crt sin 26j + d r w sin ft r + b r + c r ) 



r=l 



. K 

t sm c ^-^ _ ~ 
+ cos— — ya s d s e 

K ° 1=1. 



-*ry: 



a:atco S 2b 3 +~a s ^cosb 3 gin(a^* sin 26 a + 5 s rosin 6 S + 6 S + c s ) 
afdi(a - l) e a i t+a ^ cose 



By (5.4)-(5.6), (5.9) and (5.10), $1 = 0, 



(5.24) 



cos c 

$2 = ( —5 — / _ crRtp sin c ) cos c, 

itn 



(5.25) 
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COS c 

$3 =(——/- aRi) sine ) sine - oRz. (5.26) 
Rn 



According to (3.4), 

R cos 2 c 2t sin c 



V = 



COS — -— V rf. e «?*cos2b l +a iro cos6 l ^ 2ft . + ^ ^ ft . + ^ 



sin c i?o 



i=i 



Ro ~ 



k 

cose . tsinc t r.2 



n sin—— J2d s e d2jcos2ba+a ^ cosb3 sm(a 2 tsm2b s + h s zasmb s + c s ) 



s=l 

m 



+Rsmc J2d i e a l tcos2bi+ai ™ cosbi sin(a 2 £ sin 26* + cwsint^ + a) + ^-z 2 (5.27) 

i=l 

modulo if a — 1, and 

. n 

p = -^COS^^ V J r e«^*cos2S, + a^cosS r sin (a2^ gin 2 £. + ^ sin 6 r + c r ) 
aR Ro 

. k 

_cosc ^ _smc y^^ e aj < rtoo S 2B.+a. W co B 6. s i n (a2^ sin2 ^ + 5flC7sin & s + g s ) 

s=L 

EdiRe a * t+aiZU sin 2c cose „ , „ 2 , , „ „ R 9 . . 

— rf— ^ =— ■ + R sin c > c/;e a * t+a ' m + -z 2 , 5.28 

i=1 2i? 2 (aKl-cr)2 + J R 2 sin 2 c) ^ 2 

modulo the transformation in (1.14)- (1.1 6). 
In summary, we have: 

Theorem 5.1. Let ai,bi,Ci,a r ,b r ,c r ,d r ,a s ,b s ,c s ,d s ,c be arbitrary real numbers. De- 
note w = ycosx + z sine. We have the following solutions of the three-dimensional 
stratified rotating Boussinesq equations (1.3)-(1.7): 

u — f, v — (ft sine, w = — (p cos c, T — + z, (5.29) 

where (1) f is given in (5.18), <j) is given in (5.19), ip is given in (5.15) and p is given in 
(5.27) if a = 1; (2) f is given in (5.23), (ft is given in (5.24), i> ^ s given in (5.20) and p 
is given in (5.28) when cr ^ 1. 

Remark 5.2. By Fourier expansion, we can use the above solution to obtain the 
one depending on three arbitrary piecewise continuous functions of w. Applying the 
transformation Ti in (1.12)-(1.13) to the above solution, we get a solution involving all 
the variables t, x, y, z. 
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